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Abstract 

We use a recently improved density-matrix expansion to calculate the nuclear energy 
density functional in the framework of in-medium chiral perturbation theory. Our cal- 
, culation treats systematically the effects from l-7r-exchange, iterated l7r-exchange, and 

irreducible 27r-exchange with intermediate A-isobar excitations, including Pauli-blocking 
■ corrections up to three-loop order. We find that the effective nucleon mass M*(p) entering 

the energy density functional is identical to the one of Fermi-liquid theory when employing 
D | the improved density-matrix expansion. The strength i*v(p) °f the (Vp) 2 surface-term as 

Q ' provided by the pion-exchange dynamics is in good agreement with that of phenomenolog- 

\q ■ ical Skyrme forces in the density region po/2 < p < pq. The spin-orbit coupling strength 

F so (p) receives contributions from iterated l7r-exchange (of the "wrong sign") and from 
three-nucleon interactions mediated by 2-7r-exchange with virtual A-excitation (of the "cor- 
rect sign"). In the region around po/2 ~ 0.08 fm -3 where the spin-orbit interaction in 
nuclei gains most of its weight these two components tend to cancel, thus leaving all room 
for the short-range spin-orbit interaction. The strength function Fj{p) multiplying the 
square of the spin-orbit density comes out much larger than in phenomenological Skyrme 



forces and it has a pronounced density dependence. 
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'. 1 Introduction 

> ■ 

^ ■ The nuclear energy density functional approach is the many-body method of choice in order to 
calculate the properties of medium-mass and heavy nuclei in a systematic manner [Tj. In this 
context non-relativistic Skyrme forces [21 SI El E] have gained much popularity because of their 
analytical simplicity and their ability to reproduce nuclear properties over the whole periodic 
table within the self-consistent Hartree-Fock approximation. Another widely and successfully 
used approach to nuclear structure calculations are relativistic mean-field models [HE]. In these 
models the nucleus is described as a collection of independent Dirac quasi-particles moving in 
self-consistently generated scalar and vector mean-fields. The footprints of relativity become 
visible through the large nuclear spin-orbit interaction which emerges in that framework from 
the interplay of the strong scalar and vector mean-fields. These counteract in producing the 
(attractive) central potential but act coherently to generate the strong spin-orbit potential. 



1 Work supported in part by BMBF, GSI and the DFG cluster of excellence: Origin and Structure of the 
Universe. 
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Constraints from chiral (pion-nucleon) dynamics and the symmetry breaking pattern of QCD 
at low energies have been implemented into a pertinent relativistic point-coupling Lagrangian 
in ref. [5]. In such a more constrained scheme the results for nuclear observables come out 
comparable to those of the best phenomenological parameterizations of relativistic mean-field 
models. 

A complementary approach in the quest for predictive nuclear energy density functionals 
[TOT [TTJ [12] focusses less on the fitting of experimental data, but attempts to constrain the 
analytical form of the functional and the values of its couplings from many-body perturbation 
theory and the underlying two- and three-nucleon interactions. Switching from the conventional 
hard-core NN-potentials to low-momentum interactions is essential in this respect, because 
the nuclear many-body problem formulated in terms of low-momentum interactions becomes 
significantly more perturbative. Indeed, second-order perturbative calculations provide already 
a good account of the bulk correlations in infinite nuclear matter [13] and in doubly-magic 
nuclei [H]. 

In many-body perturbation theory the contributions to the energy are written in terms of 
density-matrices and propagators convoluted with the finite-range interaction vertices, and are 
therefore highly non-local in both space and time. In order to make such functionals numerically 
tractable in heavy open-shell nuclei it is desirable to develop simplified approximations to 
these functionals expressed in terms of local densities and currents only. At this stage of 
the procedure the density-matrix expansion comes prominently into play as it removes the 
non-local character of the exchange (Fock) contribution to the energy by mapping it onto a 
generalized Skyrme functional with density dependent couplings. Until recently, the prototype 
for that has been the density-matrix expansion of Negele and Vautherin [15] . This version of 
the density-matrix expansion (in particular its Fourier-transform to momentum space) has also 
been used in refs. [TUl ITT] to calculate the nuclear energy density functional in the framework of 
in-medium chiral perturbation theory. The density-dependent coupling strengths of the surface 
term, (Vp) 2 , or the spin-orbit term, Vp- J, arise in these calculations exclusively from the long- 
range Itt- and 27r-exchange dynamics in an inhomogeneous many-nucleon system characterized 
by a local density p(r ) and a local spin-orbit density J{r). 

In a recent paper by Gebremariam, Duguet and Bogner [IB] an improved density-matrix 
expansion has been developed for spin-unsaturated nuclei. It has been demonstrated that 
phase-space averaging techniques allow for a consistent expansion of both the spin-independent 
(scalar) part as well as the spin-dependent (vector) part of the density-matrix. A further key 
feature of the new method has been to take into account the deformation displayed by the local 
density distribution at the surface of most nuclei. The accuracy of the new phase-space averaged 
density-matrix expansion and the original one of Negele and Vautherin has been gauged via 
the Fock energy (densities) arising from (schematic finite-range) central, tensor and spin-orbit 
interactions for a large set of semi-magic nuclei. For a central force the Fock energy depends 
primarily on the spin-independent (scalar) part of the density-matrix and a few percent accuracy 
is reached for both variants of the density-matrix expansion. On the other hand the Fock energy 
due to a tensor force is determined by the spin-dependent (vector) part of the density-matrix. 
In that case the original density-matrix expansion of Negele and Vautherin leads to an error of 
about 50%, whereas the new one based on phase-space averaging techniques reduces the error 
drastically to only a few percent. This is the same level of accuracy as obtained for interaction 
terms involving the spin-independent (scalar) part of the density-matrix. For further details 
on these extensive and instructive test studies we refer to ref. [IB] • 

The purpose of the present work is to match with these new developments and to reconsider 
the nuclear energy density functional as it emerges from chiral pion-nucleon dynamics on the 
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basis of the improved density-matrix expansion of ref.[18j. Our paper is organized as follows. In 
section 2 we recall the explicit form of the improved density-matrix expansion of Gebremariam, 
Duguet and Bogner [TS]. Its Fourier-transform to momentum space provides the adequate 
technical tool to calculate the nuclear energy density functional in a diagrammatic framework. 
As a first interesting result we find that for the zero-range Skyrme force the new and the old 
density-matrix expansion lead to identical results. Differences between the two versions are 
therefore to be expected for the interaction contributions arising from the long-range In- and 
27r-exchange between nucleons. In section 3, we present the analytical results for the density- 
dependent strength functions F T (p), F so (p) and Fj(p) from which the nuclear energy density 
functional is composed. We restrict ourselves here to the isospin-symmetric case of equal 
proton and neutron number. These analytical expressions give individually the effects due to 
l7r-exchange, iterated l7r-exchange, and irreducible 27r-exchange with intermediate A-isobar 
excitations, including Pauli-blocking corrections up to three-loop order. Section 4 is devoted to 
a discussion of our numerical results and finally section 5 ends with a summary and concluding 
remarks. In the appendix the three-body spin-orbit coupling strength F so (p) is presented for 
an alternative description of the 27r-exchange three-nucleon interaction. 



2 Improved density-matrix expansion and energy den- 
sity functional 

The starting point for the construction of an explicit nuclear energy density functional is the 
density-matrix as given by a sum over the occupied energy eigenfunctions ^> a of the (non- 
relativistic) many-fermion system. According to Gebremariam, Duguet and Bogner [TB] the 
bilocal density-matrix can be expanded in relative and center-of-mass coordinates, a and f, as 
follows a 



t - -phi - -V 2 p 

ZKf - L ' 

3i 



3 %2 1. 
5 r -*' 4 



+^j^ji(ak f )a ■ (ax J) + ... , (1) 

where j\ (x) = (sin x—x cos x) / x 2 is the spherical Bessel function of index 1. The other quantities 
appearing on the right hand side of eq.(l) are the local nucleon density: 

2k 3 (r) 

P(r) = ^ 1 = ^Ur)^ a (r), (2) 
written here in terms of the local Fermi- momentum kf(r), the local kinetic energy density: 

r(0 = £v-*KO-v-*a(0. (3) 

a 

and the local spin-orbit density: 

J(f) = £tft(f)i* x W Q (f) . (4) 



2 We are considering for equal proton and neutron number the spherical phase-space averaged version without 
quadrupolar deformation of the local Fermi momentum distribution. It brings about already most of the 
improvements |18) . 
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For notational simplicity we have dropped their argument f in eq.(l) and will do so in the 
following. It is important to note that a pairwise filling of time-reversed orbitals a has been 
assumed in eq.(l), so that (various possible) time-reversal-odd fields do not come into play [1J. 
The main difference of this improved density-matrix expansion to the original one of Negele and 
Vautherin [15] lies in the index of the Bessel function multiplying the kinetic energy and spin- 
orbit densities in eq.(l). The Fourier-transform of the (expanded) density-matrix with respect 
to both coordinates a and f defines a " medium insertion" for the inhomogeneous many-nucleon 
system characterized by the time-reversal-even fields p(f), r(f) and J(r): 



T(p,q) = ld 3 re- i ^[e(k f -\p\) + ^[k f 8\k f -\p\)-28(k f -\p\\ 

x (r - \pk) - \v 2 p) - ^ 5(k f - \p\) a-(px J) j . (5) 

The double line in the left picture of Fig. 1 symbolizes this medium insertion together with the 
assignment of the out- and in-going nucleon momenta p ± q/2. The momentum transfer q is 
provided by the Fourier components of the inhomogeneous (matter) distributions p(f ), r(r ) and 
J(f ). As a check one verifies that the Fourier transform (1/27T 3 ) / d 3 pe~ ip ' a of the expression in 
the curly brackets in eq. (5) reproduces exactly the right hand side of the (improved) density- 
matrix expansion written in eq.(l). In comparison to the version of T(p,q) which followed from 
Negele and Vautherin's density-matrix expansion [16] the weight function of the kinetic energy 
density r(r ) has changed from 35(5p 2 — 3kj)9(kf— \p\) to 2k 3 [kf 8'(kf — \p\) — 28{kf — \p\)) and 

that of the spin-orbit density J(r ) has changed from 8{kf — \p |) — kf 8'{kj — \p |) to — 38(kf — \p |). 
For an inhomogeneous many-nucleon system this leads to a different weighting of the momentum 
dependent nucleon- nucleon interactions in the vicinity of the local Fermi-surface \p\ = kf(f), 
with appropriate consequences for the energy density functional. 

Going up to second order in spatial gradients (i.e. deviations from homogeneity) the energy 
density functional relevant for N = Z even-even nuclei reads: 



S[p,rJ] = pE(p) + 



3 _„i r 1 k 2 



T 



f 



+ F T {p) 



2M 4M 3 

+ (Vp) 2 F v (p) + Vp ■ JF so (p) + J 2 Fj(p) . (6) 

Here, E(p) is the energy per particle of isospin-symmetric nuclear matter evaluated at the 
local nucleon density p(f). The (small) correction term —k"j/4M 3 in eq.(6) stems from the 
relativistically improved kinetic energy and reflects in this way the relativistic increase of mass. 
The density-dependent functions F T (p), Fy(p), F so (p) and Fj(p) arising from two- and three- 
nucleon interactions encode new dynamical information specific for the inhomogeneous many- 
nucleon system. In particular, -Fy(p) measures the energy associated with density gradients at 
the nuclear surface and Fy (p) gives the strength of the spin-orbit coupling. 

Returning to eq.(5) one sees that F T (p) emerges via a perturbation on top of the density of 
states 8(kf — \p |). The single-particle potential in nuclear matter can be obtained in the same 
way by introducing a (three-dimensional) delta-function as the perturbation. Consequently, the 
strength function F T (p) can be expressed in terms of the momentum and density- dependent 
single-particle potential U(p,kf) as follows: 

1 dU(p, k f ) 

FAp) -2T f ^p— P = kf - (7) 
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Figure 1: Left: The double line symbolizes the medium insertion defined by eq.(5). Next are 
shown: The one-pion exchange Fock diagram and the iterated one-pion exchange Hartree and 
Fock diagrams. Their isospin factors for isospin-symmetric nuclear systems are 6, 12 and —6, 
respectively. 

In eq.(5) the term r — 3pk"j/5 is accompanied by — V 2 p/4. Performing a partial integration of 
the energy Jd 3 rS[p,r, J] one is lead to the decomposition: 

MP) = \ 9J j^ + F d {p) , (8) 

where F d (p) comprises all those contributions for which the (Vp) 2 -factor originates directly 
from the momentum dependence of the interactions in an expansion up to order q 2 . Since 
no information about the density-matrix expansion beyond its (fixed) nucleon matter piece 
9(kf — \p\) goes into the derivation of the strength function F d (p) the pertinent contributions 
from the 27r-exchange dynamics are still given in unchanged form by eqs. (12, 15, 19, 24, 28) in 
ref. [16] and eqs. (26,30) in ref . [IT] . 

As a first test case for the improved density- matrix expansion (summarized in eq.(5)) we 
have applied it to the (zero-range) Skyrme force [21 EES] and found that it gives identical results: 

FAp) (Sk) = ^(3fa + 5t 2 ) , F d (p)^ = 1(3*1 - 5t 2 ) , 
16 32 

F so (p) iSk) = \W , Fj(pfV = 1 2 ) , (9) 

for the energy density functional as the original density-matrix expansion of Negele and Vau- 
therin [15j . Obviously, for contact-interactions with their simple quadratic momentum depen- 
dence the different weighting of interaction strength in the vicinity of the Fermi-surface has 
no visible effect. A stronger influence of the actual form of the density-matrix expansion is 
therefore expected for the contributions arising from the long-range lir- and 27r-exchange. The 
pertinent analytical expressions are collected in the next section. 

3 Diagrammatic calculation 

In this section we present analytical formulas for the three density-dependent strength functions 
F T (p), F so (p) and Fj(p) as derived (via the improved density-matrix expansion [T8]) from l7r- 
exchange, iterated l7r-exchange, and irreducible 27r-exchange diagrams with intermediate A- 
isobar excitations, including Pauli-blocking corrections up to three-loop order. We give for 
each diagram only the final result omitting all technical details related to extensive algebraic 
manipulations and solving elementary integrals. Further explanations about the organization 
and performance of our diagrammatic calculation can be found in section 3 of ref.[16j. 
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3.1 One pion exchange Fock diagram with two medium insertions 



The non- vanishing contributions from the l7r-exchange Fock diagram shown in Fig. 1, including 
the relativistic l/M 2 -corrections, read: 



Fr{ P ) = 



87T/ W ) 2 «3 

3g\m\ 



(u 2 + -) ln(l + 4 M 2 )-2 M 2 



(Stt^M) 2 
Fj{ P ) = 



4w 3 

4u 2 arctan 2u - u ln(l + 4u 2 ) 



8u 4 - Au 2 + ln(l + 4m 2 ) 



(32m^) 2 M e 

where we have introduced the convenient dimensionless variable u = kf/m n . 



(10) 
(11) 



3.2 Iterated one-pion exchange Hartree diagram with two medium 
insertions 

The two-body contributions from the iterated l7r-exchange Hartree diagram in Fig. 1 read: 

, . q\Mml (7 + 30m 2 9x 14 1 . . 

FAp) = jkwA~^^ Hl+ ] 16arctan2M | ' ^ 



F -so(p) = , 4 Uu arctan 2u - 3u 2 - - ln(l + 4m 2 ) } . (13) 

7rm 7r (4/ 7r M) 4 I 4 J 

The expression for F so (p) in eq.(13) gives (part of) the "wrong-sign" spin-orbit interaction in- 
duced by the pion-exchange tensor force in second order. It is completed by the Fock (exchange) 
contribution and the respective Pauli-blocking corrections (see eqs.(15, 18, 22,25)). 



3.3 Iterated one-pion exchange Fock diagram with two medium in- 
sertions 

We find the following contributions from the right diagram in Fig. 1 with two medium insertions 
on non-neighboring nucleon propagators: 



Fr{ P ) = 



2x 2 - u 2 



+ 2x 2 



x 



'1 + 8x 2 + 8a; 4 ) arctanrr - (1 + 4a; 2 ) arctan 2x 



(14) 



Fso(p) 



3g\M 



27rm 7r (4^M) 4 



u 2 + 



[ dx — r 

Jo 1 + 2a; 2 



x 



4x 2 (l + x 2 ) arctanx — (1 + 4a; 2 ) arctan 2a; 



(15) 



Fj(p) = 



9g\M 



,, — / dx (u 2 - x 2 )(l + 4a; 2 ) arctan 2x 

nm^^u) 4 U 2 Jo l + 2x 2 L v A ; 

+2(m 2 - 2a; 2 + 2u 2 x 2 - 6a; 4 + 2u 2 x 4 - 6a; 6 ) arctan a;] - m 2 | . 



(16) 
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3.4 Iterated one-pion exchange Hartree diagram with three medium 
insertions 

In our way of organizing the many-body calculation, the Pauli-blo eking corrections are repre- 
sented by diagrams with three medium insertions. The corresponding contributions from the 
iterated l7r-exchange Hartree diagram read: 

3giMm 2 ( 21/ „ 2x 2m 2 r 1 , 2 , l + y 



X 



4„,2 



(i + W) 2 (6 " V + 2/2 " 2 " 2) " 4mV + ln(1 + 4 " V) 

■ 2 r 1 2uxy 



ru x * ,1 

+ dx — dy 
Jo u 6 J-i 



2s 2 + s A 
l + s 2 



It 2 — X 2 !/ 2 

(w 2 — x 2 y 2 )s 4 



+ ln 



21n(l + s 2 ) 



M 2 (l + s 



2\3 



u — xy 

((5 + s 2 )s' 2 + (s + s 3 )(s" -2s')) 



(17) 



with the auxiliary function s = + — a; 2 + x 2 ?/ 2 ^ an d its partial derivatives s' = uds/du 
and s" = u 2 d 2 s/du 2 . 



Fso{p) = 



3g\M 



K 2 m«{AUy Jo ifi 



U rj£ 

dx 



2 A 



u + xy u(5x 2 y 2 — 3w 2v 
4xy In h 



u 2 — X 2 |/2 



x 



5s + 



— 6 arctan s 



(l + s 2 ) 2 
2s 4 s'(s — 2xy) j u + xy 



u — xy 
us b (u 2 + x 2 y 2 ) 
(1 + s 2 ) 2 (u 2 - x 2 y 2 ) 



u — xy 



Fj(p) = 



(l + s 2 ) 2 
9g\Mu 3 

16n 2 m n f£ Jo ~* (1 + Au 2 y 2 ) 



dy- 



2y + (l- y 2 ) In 



2n ,_ 1 + V 



(18) 
(19) 



3.5 Iterated one-pion exchange Fock diagram with three medium 
insertions 

The evaluation of this diagram is most tedious. It is advisable to split the contributions to the 
strength functions F T (p), F so (p) and Fj(p) into "factorizable" and "non-factorizable" parts. 
These two pieces are distinguished by the feature of whether the nucleon propagator in the 
denominator can be canceled or not by terms from the product of 7riV-interaction vertices in 
the numerator. We find the following "factorizable" contributions: 



Fr( P ) 



Zg\Mml\l M A 2 n 4xi M A 2 , l + 6w 2 + 8w 4 2 2 . 

' -(1 + 4w + 2w 4 )ln(l + 4w) — — ln 2 (l + 4w) 



(4tt/^) 4 [8 
■u 4 u 2 



64m 2 



y-T + 2 ldx 



| J dx u(l + u 2 + x 2 ) - (l + (u + x) 2 ) (l + (it - x) 2 )l 



x 



u 



u 



(1 - u - x )L + u - — 

; ' l + (u + x) 2 l + (u-x) 2 

with the auxiliary function: 

1 l + (u + x) 2 
L(x, u) — — In — . 

v ' Ax l + (u-x) 2 



(20) 



(21) 
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F So{p) = 9 3g f^ si UiH 1 + 4 ^ 2 ) ~ 7U ^ arCtaI1 2U + 28 ^ 3 + 8U + ~ 

tHm^^/Vm) 4 I M 



3 + 14m 2 + 10m / 



2m 3 

+4 / daJi/ 



In 1 + 4m 2 + — — In 2 1 + 4m 2 

16m 



3a;- 2 (l + m ) + 3 + 2m 2 - m 4 - (3 + 7m 2 ):t 2 + 5x 



I 

Jo 

-6m:t- 2 (1 + m 2 ) 2 L + 3m 2 x" 2 (1 + m 2 )|| , 



(22) 



Fj(p) 



9g\M 



7r 2 m ff (8/ ff «) 4 
37m 8 - 7m 2 - 12m 4 



> arch' i! 2// lu 3 -- - 1ii 2 (1 + 4m 2 ) 



u 



16m 5 



+ 



16m 3 



1ii(1 + 4m 2 ) + J^ dx { l 



2\4 



2x 2 



;i + m 2 ) 



+2(1 - m 4 )(1 + m ) + (5 + 2m 2 + 5m 4 )x 2 - (6 + 10m 2 ):t 4 + 



lis 6 



+ - 



u 



3m 4 + 2m 2 



2\3 



ar 



+ 



2x 2 



;i + m 2 ) 2 }}. 



(23) 



The " non-factorizable" contributions (stemming from nine-dimensional principal value integrals 
over the product of three Fermi-spheres of radius kf) read on the other hand: 

3^Mm 2 /-i /-i j/ ^( j/ 2 + ^ 2 -l) f 4MV(2^ 2 -l) 
AP) ~ (47TA) 4 JJ V J-i\yz\y/y* + z*-l{ 1 + Au 2 z 2 



X 



X 



ln(l + 4MV)-4MV r _^, yo ^ 2M5(i + s2)2 
t 2 - ln(l + t 2 )} Us + s 3 )(2s' - s") - (3 + s 2 )s' 21 



9(y)0(z)+ [ U dx 
Jo 



x z s z 



(24) 



FUP) = n 2 m«{AUylj y U Z 



i yz6(y 2 + z 2 -1) (8y 2 z6( y y)6(z) 



\yz\Vy 2 



1 + 4m V 



X 



arctan(2M^) — 2uz 



ru 
/ ^ 

Jo 



+ / ofe 



x 2 sVt¥ 



2M 8 (l + S 2 )(l+t 2 ) 



(txy — srrz — si) > , (25) 



pi 9fl 4 M ri fi^ yze(y* + z*-l) f 2y*6(y)e(z) 

AP) 7r 2 m ff (4/ w ) 4 7-1 y 7-i |yz|Vy 2 + ^ 2 -l Im(1 + 4mV) 



x 



ln(l + 4mV) - 4mV 



», x 4 sWt' (1 - y 2 - z 2 ) 

rfr 

o 4M 10 (l + S 2 )(l+t 2 ) 



(26) 



with the auxiliary function t = xz + y/u 2 — x 2 + x 2 2; 2 and its partial derivative £' = udt/du. 
For the numerical evaluation of the dy dz-double integrals in eqs. (24,25,26) it is advantageous 
to first antisymmetrize the integrands in y and z and then to substitute z = \/y 2 ( 2 + 1 — y 2 - 
This way the integration region becomes equal to the unit-square < y, £ < 1. 



3.6 Irreducible two-pion exchange 

At next order in the small momentum expansion comes the irreducible 27r-exchange including 
(also) intermediate A-isobar excitations. We employ a (subtracted) spectral-function represen- 
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tation of the 7riVA-loops and find the following (two-body) contributions: 



f Ap) = 7T^ / dfi lm(V c + 3W C + 2p 2 V T + 6p 2 W T ) 

QTT' 3 J2m 7T 



2p Ztf_J L( , e) 
k f + 3p? 2B ^ + f) 



Ak 2 

In ( 1 + —f 

p 2 



(27) 



Fj(p) 



+lm(V T + 3W T ) 



dp I lm(V c + 3W ( 



c) 



1L 

Z..2 



P 

Ak) 

± ji_ 

3/i + 2k) 



(ft 2 + 2k 2 ) In ( 1 + 



4^ 
P 2 



8k 



JL 
k) 



ii 3 / Ak 2 

" (p 2 + 4k 2 f ) In ( 1 + ^ 



(2? 



The imaginary parts ImVc, ImWc, ImVr and ImPVr of the isoscalar and isovector central and 
tensor NN-amplitudes due to 27r-exchange with single and double A-excitation can be found in 
section 3 of ref.[20j. The additional contributions from the irreducible 27r-exchange with only 
nucleon intermediate states are accounted for by inserting into eqs. (27,28) the imaginary parts: 



ImW, 



c 



Am 2 . 



3^(4^)4 



Aml{l + Ag 2 A 



&9a) 



l i 2 {23g\ - Wg\ - 1) + 



48^>t 
Ami 



(29) 



ImVr 



q 9a\ P 2 -4m2 



(30) 



vr/i(4/.) 4 

At leading order the irreducible 27r-exchange generates no spin-orbit NN-interaction. It emerges 
first as a relativistic 1/M-correction. In order to see the size of such relativistic effects we have 
evaluated the energy density functional with a two-body interaction composed of the (isoscalar 
and isovector) spin-orbit NN-amplitudes Vso and Wso written in eqs. (22,23) of ref.[2T]. We 
find with it the following contribution to the spin-orbit coupling strength: 



Fso(p) 



nM(AU) 



9a- s 



- ln(l + u 2 ) - - 



18 3g 2 A , (2g\ 



+ 



4 



\2u\ 



u 



arctan u 



(31) 



which has been subtracted at p = in order to eliminate (regularization dependent) short- 
distance components. As a consequence of that subtraction only the Fock terms are included 
in the expressions in eqs. (27,28,31). 



3.7 Three-body diagrams with A-excitation 

The Pauli-blocking correction to the 27r-exchange with single A-excitation is equivalent to the 
contribution of a (genuine) three-nucleon force. In fact, one is dealing here with the same three- 
nucleon interaction as originally introduced by Fujita and Miyazawa [22J. Moreover, it has been 
shown in ref.[T3] that the inclusion of this long-range 27r-exchange three-nucleon interaction is 
essential in order to reproduce the empirical saturation point of nuclear matter when using 
the low-momentum NN-potential V[ ow _k in Hartree-Fock calculations. It is therefore equally 
interesting to see its effects on the nuclear energy functional £[p,r, J]. 

The pertinent Hartree and Fock three-body diagrams related to 27r-exchange with virtual 
A-excitation are shown in Fig. 2. The central diagram with parallel pion-lines vanishes for 
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Figure 2: Hartree and Fock three-body diagrams related to 27r-exchange with virtual A-isobar 
excitation. Their isospin factors for isospin-symmetric nuclear systems are 8, 0, and 8, respec- 
tively. 

isospin-symmetric nuclear systems. Returning to the medium insertion written in eq.(5) we 
find from the left three-body Hartree diagram in Fig. 2 the following contributions: 
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(32) 

(33) 
(34) 



with A = 293 MeV the delta-nucleon mass splitting. We have used the value 3/ a/2 for the 
ratio between the tcNA- and 7riViV-coupling constants. Note that the expression for F so (p) 
in eq.(33) gives the (dominant part of the) three-body spin-orbit coupling strength suggested 
originally by Fujita and Miyazawa [22J. The three-body effects on the energy density functional 
are completed by the contributions from the right Fock diagram in Fig. 2 which read: 
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with the auxiliary functions: 

G s (x,u) 
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(36) 
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In the appendix we present the three-body spin-orbit coupling strength F so (p) for an alternative 
description of the 27r-exchange three-nucleon interaction (using 7r7riV./V-contact vertices instead 
of propagating A-isobars). A good check of all formulas collected in this section is provided 
by their Taylor series expansion in kf. Despite the superficial opposite appearance the leading 
term in the /^-expansion is always a non- negative power of kf (which is higher for three-body 
contributions than for two-body contributions). 

4 Results and discussion 



In this section we present and discuss our numerical results obtained by summing the series of 
contributions given in section 3. The physical input parameters are: qa = 1-3 (nucleon axial 
vector coupling constant), f n = 92.4 MeV (pion decay constant), m n = 135 MeV (neutral pion 
mass) and M = 939 MeV (nucleon mass). We recall that with these physical parameters and 
a few adjustable short-distance couplings the nuclear matter equation of state E(p) and many 
other nuclear matter properties [T7j can be well described by the chiral pion-nucleon dynamics 
treated to three-loop order. 
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Figure 3: The effective nucleon mass M*(p) divided by the free nucleon mass M as a function 
of the nuclear density p. 

Returning to the energy density functional S[p,r,J] in eq.(6) one observes that the ex- 
pression multiplying the kinetic energy density r(r ) has the meaning of a reciprocal density- 
dependent effective nucleon mass: 

(40) 

It is identical to the so-called " Landau" mass introduced in Fermi-liquid theory, since it derives 
in the same way from the slope of the single-particle potential U(p,kf) at the Fermi-surface 
p = kf. This consistency of effective nucleon masses follows from the improved density- matrix 
expansion of Gebremariam, Duguet and Bogner [18], but it did not hold for the original density- 
matrix expansion of Negele and Vautherin [T5J [16] . 

Fig. 3 shows the ratio of effective to free nucleon mass M*(p)/M as a function of the nuclear 
density p = 2k^/3ir 2 . One observes a reduced effective nucleon mass which reaches the value 
M*(po) = 0.874M at nuclear matter saturation density po — 0.16 fm~ 3 . This is compatible 
with the range 0.7 < M*(p )/M < 1 spanned by phenomenological Skyrme forces [31 IH El E]- 
Somewhat unusual is the non-monotonic progression of the curve in Fig. 3. It reveals that 
a sufficiently high density (0.4p ) has to be reached until the subleading 7riVA-dynamics can 
revert the tendency of the iterated l7r-exchange to increase the effective nucleon mass. The 
same feature has also been observed for the p-wave Landau parameter f\(kf) in ref. [23] (see 
Fig. 3 therein), a quantity which is intimately related to the effective nucleon mass M*(p). 

Next, we show in Fig. 4 the strength function -Fv(p) belonging to the (Vp) 2 surface-term. 
The dashed line corresponds to the truncation to In- and iterated l7r-exchange, whereas the 
full line includes in addition the 27r-exchange and the associated three-body contributions. 
Taking the band spanned by phenomenological Skyrme forces [31 IH El [6] as a benchmark one 
may conclude that the subleading 27r-exchange dynamics leads to some improvement. The 
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Figure 4: The strength function -Fv(p) of the surface-term (Vp) 2 in the nuclear energy density 
functional versus the nuclear density p. The dashed line corresponds to the truncation to l7r- 
and iterated l7r-exchange. The full line includes also 27r-exchange and associated three-body 
contributions. 

improved density-matrix expansion [TS] has furthermore flattened and shifted downward the 
curve for Fy(p) in comparison to our previous calculation (see Fig. 8 in ref.[T7j) based on the 
Negele-Vautherin density-matrix expansion. We also note that in the relevant density region 
Po/2 < p < po the main contribution to the strength function Fy(p) comes from the component 
Fd{p) (see eq.(8)) which is insensitive to the density-matrix expansion beyond its fixed nuclear 
matter part 9(kf — \p\). 

Of particular interest is the strength F so (p) of the spin-orbit coupling provided by the explicit 
pion-exchange dynamics. The dashed curve in Fig. 5 shows the "wrong-sign" spin-orbit coupling 
strength arising from iterated l7r-exchange (i.e. the pion-exchange tensor force in second order). 
Its value at half nuclear matter density p /2 = 0.08 fm -3 decomposes as ((—124.4 + 47.7) + 
(76.5 — 35.0)) MeVfm 5 = — 35.1 MeVfm 5 into Hartree and Fock pieces supplemented by the 
respective Pauli-blo eking corrections. This net negative result amounts to about —40% of 
the empirical spin-orbit coupling strength F^ mp ^ ~ 90 MeVfm 5 . In comparison to our previous 
calculation [16J based on the Negele-Vautherin density-matrix expansion which gave at po/2 the 
value -58.1 MeVfm 5 the " wrong-sign" spin-orbit coupling strength has substantially decreased 
in magnitude. The full line in Fig. 5 shows the spin-orbit coupling strength after including the 
subleading 27r-exchange, in particular the three-body contributions eqs. (33,38). One finds now 
a pronounced cancellation in the density region around po/2 = 0.08 fm~ 3 , where the spin-orbit 
interaction in nuclei gains actually most of its weight. Such an almost complete cancellation 
leaves then all room for the short-distance NN-dynamics (not treated explicitly in this work) 
to account for the strong spin-orbit coupling in nuclei. In fact, it has been shown in ref. [24J 
that the empirical value F^ mp ^ ~ 90 MeVfm 5 of the spin-orbit coupling strength in nuclei is 
in perfect agreement with the one extracted from realistic nucleon-nucleon potentials. The 
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Figure 5: The strength function F so (p) of the spin-orbit coupling term Vp • J in the nuclear 
energy density functional versus the nuclear density p. Dashed line: In- and iterated lir- 
exchange only. Full line: 27r-exchange and three-body contributions added. 

intimate connection between the strong Lorentz scalar and vector mean-fields and the (short- 
range) spin-orbit part of the NN-potential has been elucidated in ref . [25j via (relativistic) Dirac- 
Brueckner calculations of the in-medium nucleon self-energy. 

Moreover, we note that the spin-orbit coupling strength generated by the irreducible 2tt- 
exchange as a relativistic 1/M-correction (see eq.(31) in section 3.6) contributes little to the can- 
cellation between "wrong-sign" and "correct-sign" parts shown in Fig. 5. At po/2 = 0.08 fm -3 
this piece amounts to just about —4.0 MeVfm 5 . Furthermore, we have convinced ourselves that 
the spin-orbit NN-amplitudes from 27r-exchange with A-excitation (Vso and W$o collected 
in the appendix of ref.[20j) lead to an even smaller effect. These NN-amplitudes make up a 
two-body contribution to F so (p) that scales again with 1/M. 

Finally, we show in Fig. 6 the strength function Fj(p) belonging the squared spin-orbit 
density J 2 in the nuclear energy density functional as a function of the nuclear density p. One 
observes that the inclusion of the sub leading 27r-exchange strongly reduces the values of Fj(p). 
In comparison to the (narrow) band spanned by phenomenological Skyrme forces [31 ffl El E] our 
prediction for the strength function Fj(p) is much larger in the whole density region < p < po- 
In addition, the density dependence of Fj(p) comes out markedly different, due to the long-range 
character of the pion-exchange interactions. For orientation, we reproduce by the dashed-dotted 
line in Fig. 6 the leading contribution from the l7r-exchange Fock diagram (see eq.(ll)). We 
also note that in comparison to the calculation based on the Negele-Vautherin density-matrix 
expansion the magnitude of the strength function Fj(p) has substantially increased (see Fig. 5 
in ref.|16j). 

Besides representing the non-local Fock contributions from tensor forces etc. in the energy 
density functional the J 2 -term leads to another interesting side effect. Namely, it gives rise 
to an extra spin-orbit single-particle mean-field 2Fj(p) J in addition to the "normal" one, 
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Figure 6: The strength function Fj(p) accompanying the squared spin-orbit density J 2 in the 
nuclear energy density functional versus the nuclear density p. Dashed line: Itt- and iterated 
l7r-exchange only. Full line: 27r-exchange and three-body contributions added. 

— * 

F so {p) Vp. It would be interesting to investigate the role of this additional (nucleus-dependent) 
spin-orbit mean- field together with the large values and the strong density dependence of Fj(p) 
as predicted by in-medium chiral perturbation theory. 

5 Summary and concluding remarks 

In this work we have used the recently improved density-matrix expansion of Gebremariam, 
Duguet and Bogner [T5] to calculate the nuclear energy density functional £[p, r, J] relevant for 
N = Z even-even nuclei in the framework of in-medium chiral perturbation theory. Our calcu- 
lation treats systematically the effects from l7r-exchange, iterated l7r-exchange, and irreducible 
27r-exchange with intermediate A-isobar excitations, including Pauli-blocking corrections up to 
three-loop order. 

We find that the effective nucleon mass M*(p) entering the energy density functional be- 
comes identical to the one of Fermi-liquid theory when employing the improved density-matrix 
expansion. The strength i*v(p) of the (Vp) 2 surface-term as provided by the pion-exchange 
dynamics is in good agreement with that of phenomenological Skyrme forces in the density 
region p /2 < p < p . 

The spin-orbit coupling strength F so (p) receives contributions from iterated l7r-exchange (of 
the "wrong sign") and from three- nucleon interactions mediated by 27r-exchange with virtual 
A-excitation (of the "correct sign"). In the region around po/2 — 0.08 fm~ 3 where the spin-orbit 
interaction in nuclei gains most of its weight these two components tend to cancel, thus leaving 
all room for the short-range spin-orbit interaction. The empirical value F^ mp ^ ~ 90 MeVfm 5 of 
the spin-orbit coupling strength in nuclei agrees perfectly with the one extracted from the short- 
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Figure 7: Three-body spin-orbit coupling strength F so (p) as a function of the nuclear density 
P- 

range spin-orbit component of realistic NN-potentials [21]. This part of the NN-interaction 
drives at the same time the strong Lorentz scalar and vector mean-fields on which the whole 
success of the relativistic Dirac phenomenology rests. 

The strength function Fj(p) multiplying the squared spin-orbit density comes out much 
larger than from phenomenological Skyrme forces and it has a pronounced density dependence 
due to the long-range character of the pion-exchange interaction. The interplay between the two 
components of the total nuclear spin-orbit mean-field 2Fj(p) J + F so (p) Vp should be further 
explored together with the large values and strong density dependence of Fj(p) as predicted 
by in-medium chiral perturbation theory. 

In comparison to refs. [T6"j [T7] where the density-matrix expansion of Negele and Vautherin 
has been employed, we find an improved description of the nuclear energy density functional 
S[p,T,J] on the basis of the improved density- matrix expansion [18j. In view of the fact 
that short-range contributions do not change (as exemplified here for the Skyrme force), a 
cancellation of the net two-pion exchange spin-orbit coupling strength around half nuclear 
matter density po/2 = 0.08 fm -3 is more satisfactory than having this cancellation around 
Po as discussed in ref.[2S]- In any case, the effective field theory formulation of nuclear forces 
provides short-range contributions to all four strength functions F T (p), Fy(p), F so (p) and Fj(p) 
and these can be fine-tuned in nuclear structure calculations. 

Appendix: Three-body spin-orbit coupling strength 

In this appendix we present and discuss the result for the three-body spin-orbit coupling 
strength F so (p) one obtains from an alternative description of the 27r-exchange three- nucleon 
interaction. Instead of the sequential ttN —>■ A — > nN transition with intermediate A-isobar 
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excitation one can employ the second order chiral 7T7riViV-contact vertex 



■^[25 ab (c 3 q a ■ q b - 2c x m\) + c 4 e abc T c a ■ (q a x $,)} 

J TT 



(41) 



to built up the 27r-exchange three-nucleon interaction. Here, q a ^ denote out-going pion momenta 
and we have already dropped the C2 term proportional to the product of two pion energies. In 
the present application these (off-shell) pion energies are equal to differences of nucleon kinetic 
energies, thus producing a relativistic l/M 2 -correction. The pertinent in-medium diagrams are 
those shown in Fig. 2 with the A-propagator shrunk to a point. Returning to the medium 
insertion written in eq.(5) we find from the corresponding three-body Hartree diagram (with 
two closed nucleon rings) the following contribution to the spin-orbit coupling strength: 
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where u = kf/m n . It is completed by the contribution of the three-body Fock diagram (with a 
single closed nucleon ring) which reads: 
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Fig. 7 shows the total three-body spin-orbit coupling strength F so (p) as a function of the 



nuclear density p = 2/c 3 /37r 2 for the choice of low-energy constants: 



Cl 



-0.81 GeV- 



c 3 = — 3.2GeV _1 and C4 = 5.4 GeV -1 . Its value at half nuclear matter density, F so (p /2) = 
47.5 MeVfm 5 , is now about 14% larger than the analogous A-driven three-body effects pre- 
sented in section 3.7. This small increase comes mainly from the fact that the low-energy 
constant — c 3 = 3.2 GeV -1 is not completely saturated by its dominant A-resonance contribu- 
tion #i/2A = 2.9 GeV" 1 . 

These considerations indicate that there is some uncertainty for the density at which the 
cancellation between " wrong-sign" and " correct-sign" spin-orbit coupling strength actually hap- 
pens. Nevertheless, the general features of such a balance remain unchanged. 
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